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Abstract.  Five  Anopheles  marajoara  Galvao  and  Damasceno  populations,  representing  diverse  ecological  conditions, 
were  sampled  throughout  Colombia  and  analyzed  using  nine  hypervariable  DNA  microsatellite  loci.  The  overall  genetic 
diversity  ( H  =  0.58)  was  lower  than  that  determined  for  some  Brazilian  populations  using  the  same  markers.  The  Caqueta 
population  (Colombia)  had  the  lowest  gene  diversity  ( H  =  0.48),  and  it  was  the  only  population  at  Hardy-Weinberg 
equilibrium.  Hardy-Weinberg  disequilibrium  in  the  remaining  four  populations  was  probably  caused  by  the  Wahlund 
effect.  The  assignment  analyses  showed  two  incompletely  isolated  gene  pools  separated  by  the  Eastern  Andean  cordil¬ 
lera.  However,  other  possible  geographical  barriers  (rivers  and  other  mountains)  did  not  play  any  role  in  the  moderate 
genetic  heterogeneity  found  among  these  populations  (FST  =  0.069).  These  results  are  noteworthy,  because  this  species  is 
a  putative  malaria  vector  in  Colombia. 


INTRODUCTION 

Colombia  has  ecological  conditions  over  85%  of  its  land  area 
suitable  for  permanent  malaria  transmission.12  Anopheles  dar- 
lingi  is  the  primary  malaria  vector  in  Colombia,  but  the  occur¬ 
rence  of  a  sympatric  cryptic  species  such  as  An.  marajoara ,  Galvao 
and  Damasceno,  adds  additional  complexity  to  the  study  of  the 
dynamics  of  malaria  transmission  in  this  country.  An,  marajoara 
plays  a  local  role  in  malaria  transmission  and  can  be  defined  as  a 
secondary  vector  with  a  local  importance  in  Colombia,  although 
it  is  certainly  not  a  primary  vector  like  An.  darlingi. 

The  geographic  distribution  of  An.  marajoara  ranges  from 
Costa  Rica  (Central  America)  to  Bolivia  (South  America), 
including  Trinidad  (West  Indies)  and  the  major  part  of  the 
Brazilian  Amazon  basin  (Figure  1).  An.  marajoara  has  been 
found  infected  by  Plasmodium  falciparum?  was  experimentally 
infected  with  P.  vivax,4  and  showed  physiological  resistance  to 
organochlorine  dichlorodiphenyltrichloroethane  (DDT)  insec¬ 
ticide.5  Additionally,  An.  marajoara  was  shown  to  be  the  most 
abundant  and  widely  distributed  anopheline  in  a  lot  of  urban 
areas  in  Colombia,  with  an  apparently  high  adaptation  to  human 
environments.6  As  such,  it  could  be  playing  an  important  role 
as  secondary  malaria  vector  in  some  areas  in  Colombia  and 
Brazil.3-6-9  In  fact,  in  the  Amapa  state  in  northeastern  Brazil, 
this  was  the  primary  malaria  vector.9  Additionally,  An.  mara- 
rajoara  is  also  a  possible  vector  of  Wuchereria  bancrofti  and 
Mansonella  ozzardi.10  There  is  perhaps  a  risk  of  transmitting  W. 
bancrofti  but  in  Guyana  only,  because  the  extent  of  this  para¬ 
site  has  been  greatly  reduced  in  the  neotropical  region. 

An.  marajoara  belongs  to  the  An.  albitarsis  complex,  which 
also  includes  An.  albitarsis,  An,  deaneorum  Rosa-Freitas, 
An.  oryzalimnetes  Wilkerson  and  Motoki,  An.  janconnae 
Wilkerson  and  Sallum,  and  An.  albitarsis  F.10  An,  marajoara 
and  An.  albitarsis  F  are  the  only  ones  recorded  for  this  com¬ 
plex  in  Colombia.11 

Previously,  three  genetic  studies  have  been  carried  out  with 
An.  marajoara  in  Brazil.  The  first  used  allozymes  and  restriction 
fragment-length  polymorphism  (RFLP)  of  mtDNA.  Three  allo- 
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patric  populations  were  detected,  with  a  noteworthy  heterozy¬ 
gote  deficit  and  unique  haplotypes  in  each  of  the  populations 
studied.12  The  second  study,  using  the  whole  mtCOI  gene,13 
found  multiple  signatures  of  a  demographic  expansion  event  in 
the  eastern  Brazilian  Amazon.  More  recently,  a  third  study  was 
done.  The  population  genetic  structure  of  An.  marajoara  in  15 
populations  in  Brazil  and  1  population  in  Trinidad  was  described 
using  homologous  microsatellite  markers  (J.  Conn  and  others, 
unpublished  data).  A  high  degree  of  genetic  heterogeneity  was 
found  among  the  populations  studied  (FST  =  0.198-0.401),  which 
suggested  some  incipient  allopatric  speciation  with  partial  bar¬ 
riers  to  gene  flow.  The  observed  heterozygosity  (H  =  0.604) 
was  lower  than  the  expected  heterozygosity  (H  =  0.772),  which 
could  be  explained  by  the  Wahlund  (=  subdivision)  effect. 

In  Colombia,  only  one  previous  genetic  study  has  been  con¬ 
ducted  with  this  species.  Analysis  of  42  polymorphic  random 
amplified  polymorphism  DNA  (RAPD)  markers  showed  an 
overall  expected  heterozygosity  of  0.294,  moderate  genetic 
heterogeneity  among  populations  (<()  =  0.083),  and  gene-flow 
estimates  ranging  from  1.5  to  7.5  individuals  per  generation.14 

The  goals  of  the  present  study  were:  (1)  to  determine  the  gene 
diversity  levels  and  the  Hardy-Weinberg  equilibrium  status  by 
means  of  nine  homologous  DNA  microsatellite  markers  short 
tandem  repeat  of  polymorphism  (STRP)  of  five  Colombian  An. 
marajoara  populations  distributed  in  different  geographic  and 
ecological  areas  of  this  country,  (2)  to  determine  the  degree  of 
genetic  heterogeneity  among  these  populations  and  correlate 
this  heterogeneity  with  geographic  barriers,  such  as  the  Andean 
mountains,  large  aquatic  environments,  like  the  Magdalena  and 
the  Cauca  rivers,  the  Macarena  Mountains,  and  different  eco¬ 
logical  habitats,  and  (3)  to  determine  with  different  assignment 
population  analyses  if  the  multi-locus  genotypic  profiles  of  each 
one  of  the  mosquitoes  analyzed  were  correctly  assigned  to  its 
own  population.  This  latter  aim  is  especially  important  to  assess 
potential  interchange  of  genes  of  public-health  importance,  such 
as  those  that  confer  insecticide  resistance  or  refractoriness  to 
infection  by  human  Plasmodia. 

MATERIALS  AND  METHODS 

Samples  and  molecular  analyses.  Five  populations  distributed 
across  the  geographic  range  of  An.  marajoara  in  Colombia 
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Figure  1.  Map  with  the  overall  distribution  of  An.  marajoara  in  Latin  America. 
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Figure  2.  Map  with  the  An.  marajoara  populations  sampled  in  Colombia.  Altitude  is  in  meters. 


were  sampled  (Figure  2). These  populations  were  from  Caceres, 
Antioquia  Department,  07°34'N,  75°2TW,  100  m  above  sea 
level  ( N  =  38  individuals  analyzed);  El  Banco,  Magdalena 
Department,  09°00'N,  73°58'W,  35  m  above  sea  level  ( N  =  36); 
El  Zulia,  Norte  de  Santander  Department,  07°56'N,  72°36'W, 
200  m  above  sea  level  (N  =  40);  Villavicencio,  Meta  Depart¬ 
ment,  04°09'W,  73°38'W,  467  m  above  sea  level  (N  =  39); 
and  Cartagena  del  Chaira,  Caqueta  Department,  01°21'N, 
74°50'W,  350  m  above  sea  level  ( N  =  40).  All  localities  are 
situated  on  the  banks  of  large  rivers.  Caceres  is  placed  at 
the  lower  Cauca  near  to  the  western  area  of  the  San  Lucas 
Mountain  and  between  the  Nechi  and  Cauca  rivers  with 
a  mean  temperature  of  28°C.  The  main  human  economic 
activities  are  the  gold  extraction  and  livestock  production. 
El  Banco  is  placed  at  the  confluence  of  the  Magdalena  and 
Cesar  rivers.  The  mean  temperature  is  around  29°C,  but  in  the 
hottest  period,  the  temperature  could  be  near  to  40°C.  The 
clime  is  tropical  with  wet  and  dry  periods.  The  main  economic 
activities  are  based  on  agriculture  and  fishing.  El  Zulia  is  at 
14  km  from  the  city  of  Cucuta.This  place  has  an  average  annual 
temperature  of  28°C,  and  its  main  economic  activities  are  those 
related  to  rice  and  sugar  as  well  as  mineria  (calizas,  carbon, 
y  arena).  Villavicencio  is  the  capital  of  the  Meta  Department 
and  is  near  to  the  Guatiquia  River.  The  mean  temperature  is 
27°C,  and  the  main  activities  are  livestock  (one  of  the  most 
important  areas  for  this  activity  in  Colombia)  and  petroleum 
extraction.  Cartagena  del  Chaira  is  placed  at  the  Caguan  River 
in  the  Colombian  Amazon  (neotropical  rain  forest)  with  an 
average  annual  temperature  of  30°C.  The  main  human  activity 
is  the  agriculture  (bananas,  maize,  and  yucca)  and  fishing. 
These  sites  have  relative  humidity  >  80%. The  collections  were 
based  on  accessibility  to  the  localities  and  confirmation  of  local 
security  because  of  the  question  of  the  guerrillas  in  Colombia. 

All  specimens  ( N  =  193)  were  collected  between  July  and 
September  2004.  These  were  identified  as  An.  marajoara  using 


morphological  taxonomic  keys15  and  were  then  maintained  in 
absolute  ethanol  at  4°C.  All  the  specimens  analyzed  (except 
three  individuals  sampled  in  the  Magdalena  population  that 
were  identified  as  An.  triannulatus)  were  confirmed  as  An.  mara¬ 
joara  using  inter  space  2  ribosomal  DNA  (ITS2)  (rDNA)16  and 
the  white  gene  sequences.17  DNA  was  extracted  by  a  phenol- 
chloroform  method18  that  yielded  approximately  2  pg  of  DNA 
per  specimen.  The  STRP  markers  used  were  those  designed 
by  Li  and  others19  specifically  for  the  species  studied  here,  and 
they  worked  perfectly:  M9-2,  Ml-31 ,  M8-2-2,  M2 ,  M7-1 ,  M6-8, 
M5-71,  M10-1,  and  Mll-4.  One  of  the  primers  from  each  pair 
was  labeled  with  fluorescein  amidite  (FAM),  5'hexachloro- 
flurescein-CE  phosphoramidite  (FIEX),  or  fluorescein  benzo- 
xanthene  (NED)  (AB  Prism,  Washington,  DC).  Touch-down 
polymerase  chain  reaction  (PCR)  was  conducted  in  a  20-pL 
reaction  containing  1  x  buffer  (10mMTris-HCl,pFI  8.3,50mM 
KC1),  2.5  mM  MgCl2,  200  pL  dideoxy  nucleoside  triphos¬ 
phates  (dNTP),  5  pmol  each  primer,  0.5  U  Taq  polymerase 
(AB),  and  5  ng  DNA  template.  PCRs  were  performed  on  a 
GeneAmp  PCR  System  9700  (AB).  With  the  exception  of 
M5-71,  the  reaction  conditions  were  3  min  at  95°C  followed 
by  two  cycles  of  94°C  for  30  s,  60.7°C  for  30  s  to  annealing 
temperature,  and  72°C  for  30  s;  14  cycles  using  a  touchdown 
protocol  (94°C  for  15  s,  15  s  at  58.6°C,  temperature  decreased 
by  1.7°C  every  other  cycle  to  a  touchdown  annealing  temper¬ 
ature,  72°C  for  3  s)  were  followed  by  24  cycles  with  a  stan¬ 
dard  protocol  (94°C  for  15  s,  50°C  as  annealing  temperature 
for  15  s,  and  72°C  for  30  s)  and  were  closed  with  10  min  at 
72°C  to  reduce  stutter  bands.  The  annealing  temperature  of 
M5-71  marker  was  66.7°C  for  the  first  cycle,  66°C  for  the  sec¬ 
ond  cycle,  and  55°C  for  the  third  cycle.  Products  were  mixed 
(Mix  1:  M9-2,  Ml-31,  M8-2-2,  M2,  and  M7-7;  Mix  2:  M6-8, 
M5-71,M10-1,  and  Mll-4)  and  run  in  an  AB  13100  (AB).  Data 
were  automatically  collected  and  analyzed  by  Genotyper 
software  3.7v  (AB).  The  samples  that  showed  ghost  bands  or 
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rare  alleles  (alleles  that  were  only  found  in  one  of  the  popula¬ 
tions  studied)  were  run  a  second  time  to  confirm  the  results 
obtained. 

Population  genetics  analyses.  The  mean  number  of  alleles 
per  locus  and  the  expected  heterozygosity  ( H ) 20  were 
calculated  for  each  of  the  populations,  and  differences  among 
these  estimates  were  statistically  analyzed  with  a  Student’s 
t  test. The  heterozygosity  data  were  arcsine  transformed  before 
analysis,  as  proposed  by  Archie.21  This  analysis  will  permit  us 
to  determine  if  some  of  the  populations  analyzed  have  lost 
gene  diversity  by  gene  drift  or  founder  effect. 

The  Hardy-Weinberg  equilibrium  (HWE)  was  estimated 
using  the  Weir  and  Cokerham22  1984  F  (W-C  F)  statistic. 
This  procedure  was  used  to  calculate  the  excess  of  homo- 
and  heterozygous  genotypes  within  each  population,  and  the 
Markov  chain  method  with  1,000  dememorization,  50  batches 
at  1,000  iterations  per  batch,23  was  used  to  measure  the  prob¬ 
abilities.  Fisher’s  method  was  used  to  analyze  the  global  HWE 
estimates  (by  loci  and  populations)  with  GenePop  3.4  ver¬ 
sion.24  Genotypic  linkage  disequilibrium  was  analyzed  by 
marker  pairs  in  each  population,  and  for  each  marker  pair,  all 
of  the  populations  were  taken  together.  Markov  Monte  Carlo 
chains  with  1,000  dememorization,  50  batches,  and  1,000  itera¬ 
tions  per  batch  were  applied  to  determine  the  significance  of 
the  probabilities. 

Several  strategies  were  used  to  calculate  genetic  heteroge¬ 
neity  globally  and  among  population  pairs.  First,  the  gene  fre¬ 
quencies  of  the  nine  microsatellites  were  tested  with  Markov 
chains  with  1,000  dememorizations,  20  batches,  and  1,000  iter¬ 
ations  per  batch.  Second,  Wright’s25  1951  F  statistics,  with  Weir 
and  Cokerham’s22 1983  procedure,  were  applied.  The  standard 
deviations  of  the  F  statistics  and  the  confidence  intervals  (95% 
and  99%)  were  calculated  with  jackknifing  and  bootstrapping 
(10,000)  over  loci,  respectively.  The  significance  of  FST  (genetic 
heterogeneity)  was  calculated  with  the  G  test  (10,000  allele 
randomizations,  random  mating  assumed)  and  the  log-likeli¬ 
hood  G  test  (10,000  genotype  randomizations,  random  mat¬ 
ing  not  assumed)  for  among  populations.26  The  significance  of 
FIS  and  FIT  were  similarly  calculated  (10,000  randomizations 
of  alleles  within  and  overall  for  the  populations  analyzed). 
For  this  second  task,  the  FSTAT  2. 9. 3.2  version  software  was 
used.27 

In  a  complementary  way,  Nei’s  gene  diversity  analysis20  was 
used  to  measure  the  gene  heterogeneity  between  populations 
as  well  as  the  RST  statistic.28,29  Rho,  as  an  unbiased  estimator  of 
Rsr>  was  calculated  for  individual  loci  and  for  all  the  loci  taken 
simultaneously  for  all  the  populations  studied  and  population 
pairs  using  R  software.30  From  these  heterogeneity  statis¬ 
tics,  theoretical  gene-flow  estimates  were  obtained  by  means 
of  an  infinite  island  model31-33  and  an  ^-dimensional  island 
model.34,35  Another  gene-flow  estimate  was  obtained  with  the 
private  allele  method  of  Slatkin36  and  Barton  and  Slatkin.37 

Finally,  an  assignment  analysis  was  applied  to  the  193  multi¬ 
locus  genotype  profiles  obtained.  This  assignment  analysis  was 
applied  by  means  of  the  Structure  2.1  software.38  This  method, 
which  employs  Markov  Chain  Monte  Carlo  (MCMC)  pro¬ 
cedures  and  the  Gibbs  sampler,  uses  multilocus  genotypes 
to  infer  population  structure,  and  individuals  are  simultane¬ 
ously  assigned  to  specific  populations.  The  model  considers 
K  populations,  where  K  may  be  unknown,  and  the  individu¬ 
als  are  assigned  tentatively  to  one  population  or  jointly  to  two 
or  more  populations  if  their  genotypes  are  considered  to  be 


admixed.  The  posterior  K  probabilities  are  calculated  assum¬ 
ing  uniform  prior  values  of  K  (USEPOPINF  =  0).  The  pres¬ 
ence  of  the  most  probable  number  of  gene  pools  within  data 
considered  is  revealed  by  the  increasing  likelihood.  After 
the  most  likely  number  of  populations  is  found,  the  analysis 
is  repeated  but  with  introducing  the  model  with  prior  geo¬ 
graphic  population  information  (USEPOPINF  =  1).  Different 
parameter  conditions  were  considered.  No  admixture  model 
and  admixture  model  were  analyzed.  In  the  first  case,  the  indi¬ 
viduals  exclusively  belong  to  one  population,  whereas,  in  the 
second  model,  one  individual  could  simultaneously  belong  to 
different  populations.  In  each  case,  another  parameter  was 
analyzed:  gene  frequencies  correlated  among  the  populations 
or  independent  gene  frequencies  correlated  among  the  popu¬ 
lations.  With  correlated  gene  frequencies,  two  FST  values  were 
considered  (0.001,  the  default  value  of  the  program,  and  0.069, 
the  value  empirically  obtained  for  the  present  study).  In  the 
case  of  correlated  gene  frequencies,  two  migration  rates  (0.01 
and  0.05)  were  considered  as  well.  With  each  condition,  the 
program  was  run  with  1  million  iterations  after  a  burn  period 
of  100,000  iterations.  Each  analysis  was  performed  two  times 
with  convergent  results. 

RESUFTS 

All  mosquitoes  analyzed  belonged  to  An.  marajoara  based 
on  morphological  characters  and  ITS2  (rDNA)  and  white  gene 
sequences,  with  the  exception  of  three  individuals  sampled  in 
El  Banco,  Magdalena;  these  were  deleted  from  the  population 
genetics  analyses.  Therefore,  the  homozygote  excess  was  not 
caused  by  the  mix  of  different  albitarsis  complex  species  in  our 
sample. 

The  nine  STRP  markers  studied  were  polymorphic  globally 
and  for  each  of  the  locations  analyzed.  Globally,  M9-2  and  M5- 
71  yielded  the  highest  number  of  alleles  (10),  whereas  M8-2-2 
had  the  lowest  (6)  (Table  1).  For  populations,  the  markers  with 
the  highest  number  of  alleles  were  M5-71  (9;  Caqueta),  M5-71 
(9;  Meta),  M5-71  (7;  Norte  de  Santander),  Ml-31  and  M10-1 
(7  each;  Magdalena),  and  M5-71 ,  M6-8 ,  and  M7-1  (6  each; 
Antioquia).  In  Table  1,  the  15  rare  alleles  are  shown  (alleles 
that  were  only  found  in  one  of  the  populations  studied)  that 
were  found  in  this  study. 

The  average  number  of  alleles  per  population  was  relatively 
uniform  in  all  populations.  The  two  populations  with  values 
slightly  lower  were  Caqueta  (n  =  4.56  ±  2.07)  and  Antioquia 
(nA  =  4.89  ±  1.27);  Norte  de  Santander  (nA  =  5.33  ±  1.32)  and 
Magdalena  (nA  =  5.53  ±  1.07)  yielded  intermediate  values,  and 
Meta  (nA  =  6.11  ±  1.54)  had  the  highest  value.  Nevertheless,  no 
statistically  significant  differences  were  detected  between  any 
pair  of  populations. 

Globally,  the  expected  mean  heterozygosity  (gene  diversity) 
was  H  =  0.578,  which  is  an  intermediate  value.  Per  location,  all 
the  populations  showed  very  similar  values  (Magdalena,  H  = 
0.585;  Antioquia,  H  =  0.557;  Norte  de  Santander,  H  =  0.555; 
Meta,  H  =  0.547)  with  Caqueta  presenting  a  slightly  lower 
value  ( H  =  0.481),  which  was  significantly  lower  than  that  for 
Magdalena.  Therefore,  for  both  statistics  (mean  allele  number 
and  expected  heterozygosity),  the  Caqueta  population  had  the 
lowest  genetic  diversity.  By  marker,  M5-71  ( H  =  0.729),  M8-2-2 
(H  =  0.649),  and  M10-1  (H  =  0.622)  yielded  the  highest  levels 
of  gene  diversity,  whereas  Mll-4  ( H  =  0.371)  clearly  showed 
the  lowest  diversity. 


Table  1 

Alleles  found  in  the  total  sample  and  for  each  one  of  the  An.  marajoara  samples  analyzed  for  the  nine  STRP  used 
Locus  M9-2  M-2  Ml-31  MS-2-2  M10-1  Mil -4  MS-71  M6-8  M7-1 
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HWE.  In  Table  2,  the  HWE  results  are  shown.  Globally, 
only  one  STRP  marker  was  in  HWE  (MS-2 -2;  P  =  0.77).  The 
remaining  eight  markers  were  clearly  in  HW  disequilibrium 
by  a  significant  excess  of  homozygotes  (Weir  and  Cokerham22 
F  statistic  of  0.104-0.320;  P  =  0.0146-0.0002).  Therefore, 
the  total  Fisher  test  done  simultaneously  for  the  nine 
STRP  markers  and  the  global  sample  detected  a  very  high 
homozygote  excess  with  regard  to  the  HWE  ( % 2  =  infinite,  8  df, 
P  <  0.00000).  By  population,  only  Caqueta  was  in  overall  HWE 
(F  =  0.006;  P  =  0.219).  In  this  population,  only  one  STRP  was 
not  in  HWE  but  by  heterozygote  excess  ( M8-2-2 ),  although 
when  the  Bonferroni  correction  was  applied,  this  significance 
disappeared.  The  other  four  populations  were  globally  in  HW 
disequilibrium.  Antioquia  presented  a  homozygote  excess 
(F  =  0.065;  P  =  0.0007)  mainly  from  M9-2,  Ml-31,  and  M7-1. 
With  the  Bonferroni  correction,  only  the  last  two  markers 
presented  significant  homozygote  excess.  Magdalena  showed 
a  homozygote  excess  (F  =  0.131;  P  =  0.0004)  mainly  from  M2, 
Ml-31,  M8-2-2,  and  M7-1,  although  only  Ml-31  was  significant 
after  the  Bonferroni  correction.  Norte  de  Santander  yielded  a 
significant  homozygote  excess  (F  =  0.093;  P  =  0.0000)  mainly 
from  the  contribution  of  Ml  -31  and  M6-8.  Finally,  Meta  showed 
a  noteworthy  homozygote  excess  (F  =  0.096;  P  =  0.0000)  from 
the  considerable  contributions  of  M9-2  and  M2.  Such  as  it 
was  shown,  each  population  was  in  HW  disequilibrium  by  the 
contribution  of  different  STRP  markers. 

Linkage  disequilibrium.  Linkage-disequilibrium  analysis 
applied  to  all  the  populations  taken  together,  and  it  showed 
one  case  of  significant  disequilibrium  of  the  M9-2  and  M5-7 
pair  (x2  =  infinite,  P  =  0.0000)  after  the  Bonferroni  correction. 
The  same  analysis  by  population  showed  that  only  one 
population  (Meta)  had  two  significant  linkage  disequilibrium 
at  two  loci:  at  M9-2  and  M5-7  ( P  =  0.0000)  and  M9-2  and 
Ml -31  ( P  =  0.00212).  With  the  Bonferroni  correction,  the  latter 
disappeared.  With  this  one  exception,  we  conclude  that  the 
population  dynamics  of  all  STRP  markers  are  independent. 

Genetic  heterogeneity.  The  exact  test  with  Markov  chains 
for  genetic  heterogeneity  using  gene  frequencies  showed 
that  eight  of  nine  STRP  markers  used  presented  significant 
differences  across  the  populations  analyzed.  Only  Ml-31  did 
not  yield  significant  heterogeneity  (P  =  0.899);  meanwhile, 
the  global  set  of  nine  markers  showed  a  striking  statistically 
significant  difference  with  the  Fisher’s  test  (%2  =  infinite, 
18  df,  P  =  0.0000).  The  markers  that  showed  most  significant 
differences  between  population  pairs  by  means  of  exact  tests 
with  Markov  chains  were  M9-2  (8  cases),  Ml  1-4  (7  cases),  M5-7 
(7  cases),  and  M7-1  (7  cases),  whereas  Ml-31  (1  case)  and 
M8-2-2  (2  cases)  presented  least  significant  discriminatory 
capacity  among  the  populations  studied.  The  highest  number 
of  significant  differences  between  population  pairs  were 
Caqueta-Magdalena  (8  of  9  STRP  markers)  and  Caqueta- 
Antioquia  (7  of  9  STRP).  Therefore,  the  more  geographically 
distant  populations  were  those  with  the  highest  number  of 
significant  comparison  pairs.  Antioquia-Magdalena  (1  of  9 
STRP),  Caqueta-Norte  de  Santander  (1  of  9  STRP),  and 
Meta-Norte  Santander  (3  of  9  STRP)  were  the  populations 
that  presented  the  least  number  of  significant  comparison 
pairs.  Antioquia  and  Magdalena  are  situated  in  the  western 
area  of  the  eastern  Andean  Mountains;  meanwhile,  Caqueta, 
Meta,  and  Norte  de  Santander  are  in  the  eastern  area  of 
this  cordillera.  This  is  the  first  evidence  that  the  eastern 
Andean  Cordillera  may  play  an  important  role  in  genetically 
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Table  2 

F  values  for  each  locus  and  each  Colombian  An.  mararajoara  population  according  to  Weir  and  Cokerham22 
Caqueta  Antioquia  Magdalena  Norte  de  Santander  Meta 


Locus 

F,s 

SD 

P 

Fs 

SD 

P 

Fa 

SD 

P 

F,s 

SD 

P 

Fs 

SD 

P 

M9-2 

0.183 

0.000 

0.3442 

0.152 

0.002 

0.0146* 

0.174 

0.010 

0.1326 

0.060 

0.013 

0.1321 

0.507 

0.000 

0.0000* 

M2 

0.059 

0.013 

0.1196 

-0.182 

0.008 

0.1191 

0.132 

0.004 

0.0400* 

0.080 

0.013 

0.2137 

0.383 

0.000 

0.0000* 

Ml-31 

0.082 

0.000 

0.3944 

0.355 

0.002 

0.0077* 

0.071 

0.002 

0.0034* 

0.369 

0.000 

0.0000* 

-0.005 

0.008 

0.0517 

M8-2-2 

-0.318 

0.000 

0.0207* 

-0.069 

0.000 

1.0000 

0.250 

0.004 

0.0287* 

-0.082 

0.000 

0.9960 

-0.050 

0.008 

0.7776 

M10-1 

0.100 

0.015 

0.6317 

0.137 

0.0146 

0.3985 

0.111 

0.015 

0.2093 

-0.037 

0.000 

0.9604 

0.063 

0.017 

0.6248 

Mll-4 

0.000 

0.000 

0.0000* 

-0.100 

0.000 

0.8240 

0.093 

0.011 

0.2380 

-0.039 

0.000 

1.0000 

-0.053 

0.000 

1.0000 

M5-71 

0.011 

0.049 

0.5652 

0.071 

0.019 

0.3526 

-0.132 

0.000 

1.0000 

0.054 

0.014 

0.4971 

0.007 

0.012 

0.8547 

M6-8 

0.010 

0.000 

0.7587 

0.128 

0.021 

0.2711 

0.159 

0.000 

0.5609 

0.369 

0.003 

0.0053* 

0.108 

0.012 

0.2135 

M7-1 

0.115 

0.019 

0.4971 

0.080 

0.001 

0.0009* 

0.234 

0.002 

0.0087* 

0.042 

0.014 

0.3110 

-0.108 

0.000 

1.0000 

Total 

0.006 

0.065 

0.131 

0.093 

0.096 

P 

=  0.2199 

P  =  0.0007* 

P 

=  0.0004* 

P 

=  0.00000* 

P 

=  0.00000* 

P  =  probability;  SD  =  standard  deviation. 
*  Significant  value. 


separating,  in  a  moderate  way,  the  populations  of  An.  marajora 
in  Colombia. 

The  application  of  the  Wright’s  F  statistics  with  the  Weir  and 
Cokerham22  procedure  yielded  the  following  results  (Table  3). 
The  overall  value  of  FIT,  0.145  [confidence  interval  (Cl)  with 
10,000  bootstrapping  at  99%  =  0.064-0.221),  showed  a  sig¬ 
nificant  homozygote  excess  at  the  total  population  level. 
Individually,  all  the  loci  studied,  except  M8-2-2  ( P  =  0.77),  pre¬ 
sented  significant  homozygote  excess  at  this  population  level. 
The  average  value  of  F]S  was  slightly  lower  than  the  previous 
case,  although  it  also  showed  homozygote  excess  at  the  sub¬ 
population  level  (0.082  with  a  99%  Cl  of  0.007-0.155).  M9-2, 
Ml -31 ,  M6-8 ,  and  M7-1  showed  significant  homozygote  excess. 
These  results  corroborated  the  homozygote  excess  detected 
with  the  individual  HWE  tests.  The  genetic  heterogene¬ 
ity  measured  by  FST  was  moderately  significant  (0.069;  Cl  = 
0.025-0. 148).  This  means  that  each  population  had,  on  average, 
93.1%  of  the  genetic  diversity  found  in  the  total  Colombian 
An.  marajora  population.  All  the  loci,  except  Ml -31,  showed 
significant  heterogeneity  independently  to  consider,  or  not, 
random  breeding  (G  tests  and  log-likelihood  G  test).  The 
three  markers  with  the  highest  levels  of  genetic  heterogene¬ 
ity  were  Mll-4,  M5-7,  and  M10-1.  Nei’s20  genetic  diversity 
analysis  confirmed  similar  genetic  heterogeneity  values  as 
the  Gst  and  GST'  statistics  (0.057-0.07).  When  the  RST  statis¬ 
tic  was  used,  the  average  values  ranged  from  0.042  to  0.047, 
depending  on  the  procedures  used.  This  statistic  again  showed 


a  moderate  divergence  among  the  populations,  with  the  high¬ 
est  genetic  heterogeneity  levels  for  Mll-4  ( RST  =  0.152),  M7-1 
(Rst  =  0.117),  and  M5-7  ( RST  =  0.061).  As  in  other  analyses,  the 
most  similar  pairs  were  Antioquia-Magdalena,  Meta-Norte 
de  Santander,  and  Caqueta-Meta.  That  is,  the  most  similar 
populations  were  those  at  the  western  area  from  the  eastern 
Andean  Cordillera  and  two  pairs  of  populations  east  of  the 
mountains.  The  overall  gene-flow  estimates  from  an  infinite 
island  model  were  3.37  ( FST )  and  5.07-5.70  (RST)',  meanwhile, 
for  an  « -dimensional  island  model,  these  values  were  2.16  and 
3.24-3.65,  respectively.  The  gene-flow  estimate  from  the  pri¬ 
vate  allele  model  of  Slatkin36  was  3.6,  very  similar  to  estimates 
obtained  with  other  methods.  These  gene-flow  values  could 
be  considered  moderate,  although  the  populations  are  clearly 
related  by  gene  flow.  If  we  observe  the  gene-flow  estimates  by 
population  pair  with  the  unbiased  estimator  Rho,  they  agree 
quite  well  with  previous  analyses.  Antioquia  and  Magdalena 
are  really  one  unique  population  ( Nm  =  infinite);  meanwhile, 
Meta  and  Norte  de  Santander  were  highly  homogenized  by 
gene  flow  (Nm  =  22.86-28.79)  as  well  as  Caqueta  and  Meta 
(Nm  =  13.31-25.66)  and  Caqueta-Norte  de  Santander  (Nm  = 
10.15-10.44). The  population  pairs  most  isolated  by  lower  gene 
flow  were  Caqueta-Antioquia  (Nm  =  2.36-2.80)  and  Caqueta- 
Magdalena  (Nm  =  2.68-3.51). 

Assignment  analysis.  The  assignment  analysis  with  the 
Structure  2.1  software  showed  the  same  results  assuming 
different  hypotheses.  The  least  negative  likelihood  value 


Table  3 


Wright’s  F  statistics  values  per  locus  according  to  Weir  and  Cokerham22 


Locus 

Ft 

Ft 

PI 

P2 

F1S 

M9-2 

0.258;  P  =  0.0002* 

0.047 

0.0002* 

0.0002* 

0.221;  P  =  0.0002* 

M2 

0.104;  P  =  0.0148* 

0.022 

0.0002* 

0.0002* 

0.084;  P  =  0.0550 

Ml-31 

0.185;  P  =  0.0006* 

0.021 

0.5214 

0.2198 

0.167;  P  =  0.0016* 

M8-2-2 

-0.035;  P  =  0.7704 

0.023 

0.0002* 

0.0004* 

-0.059;  P  =  0.8914 

M10-1 

0.122;  P  =  0.0080* 

0.057 

0.0002* 

0.0002* 

0.069;  P  =  0.0962 

Mll-4 

0.320;  P  =  0.0002* 

0.330 

0.0002* 

0.0002* 

-0.015;  P  =  0.6226 

M5-71 

0.146;  P  =  0.0002* 

0.135 

0.0002* 

0.0002* 

0.014;  P  =  0.3976 

M6-8 

0.176;  P  =  0.0006* 

0.019 

0.0002* 

0.0002* 

0.160;  P  =  0.0012* 

M7-1 

0.127;  P  =  0.0044* 

0.043 

0.0002* 

0.0002* 

0.088;  P  =  0.0404* 

Total 

0.145 

0.069 

0.069 

0.069 

0.082 

Cl 

0.064-0.221 

0.025-0.148 

0.025-0.148 

0.025-0.148 

0.007-0.155 

Cl  =  99%  confidence  interval. 

*  Significant  probabilities  ( P ).  PI  =  assuming  non-random  breeding;  P2  =  assuming  random  breeding. 
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was  with  K  =  2  in  all  the  cases.  Therefore,  independent  of  the 
conditions  chosen,  two  gene  pools  seem  to  be  present  in  the 
Colombian  area  sampled.  One  gene  pool  was  composed  of 
individuals  from  Antioquia  and  Magdalena  (this  gene  pool  was 
named  Western  pool);  the  other  was  composed  of  individuals 
from  Norte  de  Santander,  Meta,  and  Caqueta  (this  gene  pool 
was  named  Eastern  pool).  These  gene  pools  are  separated  by  the 
eastern  Andean  Cordillera.  With  no  geographic  information, 
admixture  model,  correlated  frequencies  among  populations, 
F  =  0.069,  and  migration  rate  of  0.05  in  the  Western  pool,  all 
the  Antioquia  and  Magdalena  individuals  appeared,  with  the 
exception  of  three  individuals  of  Antioquia  and  five  individuals 
of  Magdalena  that  were  assigned  with  higher  probability  to 
the  Eastern  gene  pool,  plus  one  individual  from  Caqueta,  four 
individuals  from  Meta,  and  eight  individuals  from  Norte  de 
Santander  (these  last  13  specimens  sampled  in  the  Eastern 
gene  pool  but  showed  more  genetic  resemblance  with  the 
Western  gene  pool).  In  the  Eastern  gene  pool,  the  vast  majority 
of  the  individuals  of  Caqueta,  Meta,  and  Norte  de  Santander 
were  assigned,  with  the  aforementioned  exceptions  and  the 
three  and  five  individuals  from  Antioquia  and  Magdalena 
previously  mentioned.  One  hybrid  individual  between  the 
Eastern  and  Western  gene  pool  was  detected  in  Caqueta,  and 
two  hybrid  individuals  between  the  Eastern  and  Western  gene 
pool  were  found  in  Meta  and  Norte  de  Santander  (Figure  3A). 
With  other  conditions,  such  as  a  priori  geographic  information, 
admixture  model,  and  independent  gene  frequencies  among 
the  populations,  the  results  were  similar.  In  the  Eastern  gene 
pool,  all  the  individuals  from  Caqueta  and  the  major  part 
of  the  individuals  from  Meta  and  Norte  de  Santander  were 
assigned,  plus  five  individuals  from  Magdalena.  In  the  Western 
gene  pool,  all  the  Antioquia  individuals  and  the  major  fraction 
of  the  Magdalena  individuals,  together  with  four  individuals 


from  Meta  and  eight  individuals  from  Norte  de  Santander, 
were  included.  Only  one  hybrid  individual  between  the  Eastern 
and  Western  gene  pool  in  Meta  and  four  hybrid  individuals 
between  these  two  gene  pools  in  Norte  de  Santander  were 
detected  (Figure  3B).  Therefore,  the  final  picture  was  as 
follows:  two  gene  pools  (Eastern  and  Western  ones),  separated 
by  the  Eastern  Andean  Cordillera,  were  detected  but  with 
incomplete  isolation,  because  mosquitoes  from  one  gene 
pool  were  found  in  the  other  gene  pool  and  vice  versa.  The 
two  populations  that  showed  the  highest  number  of  hybrids 
between  both  Eastern  and  Western  gene  pools  were  Norte  de 
Santander  and  Meta.  Therefore,  the  two  most  differentiated 
populations  were  those  occupying  the  most  extreme  sampling 
geographic  positions:  Caqueta  representing  the  Eastern  gene 
pool  and  Antioquia  representing  the  Western  gene  pool. 

DISCUSSION 

The  overall  genetic  diversity  found  in  Colombian  An.  mara- 
joara  populations  (H = 0.578)  was  relatively  modest  with  regard 
to  that  determined  for  Brazilian  populations  of  An.  marajoara, 
H  =  0.715  (J.  Conn  and  others,  unpublished  data).  However, 
the  populations  from  Trinidad  ( H  =  0.538)  and  Boa  Vista 
( H  =  0.486),  also  studied  by  the  same  authors,  presented  values 
very  similar  to  those  reported  here.  These  values  suggest  that 
central  Brazilian  and  eastern  Amazon  populations  may  repre¬ 
sent  the  central  range  of  the  species,  whereas  the  Colombian 
populations  were  probably  founded  in  a  subsequent  expan¬ 
sion  with  a  certain  loss  of  gene  diversity  because  of  founder 
effect.  In  fact,  Lehr  and  others39  determined  that  the  An.  mara¬ 
joara  populations  underwent  a  significant  fragmentation  event 
during  the  Pleistocene.  All  the  Colombian  populations  ana¬ 
lyzed  showed  very  similar  levels  of  expected  heterozygosity 
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Figure  3.  Assignment  analysis  of  the  An.  marajoara  individuals  studied  by  means  of  two  Structure  2.1  v  software  procedures.  (A)  Assignment 
analysis  with  no  geographic  information,  admixture  model,  correlated  frequencies  among  populations,  F  =  0.069,  and  migration  rate  of  0.05.  (B) 
a  priori  geographic  information,  admixture  model,  and  independent  gene  frequencies  among  the  population.  Kl  =  western  Colombian  gene  pool; 
K2  =  eastern  Colombian  gene  pool;  H  =  hybrid  individuals  found. 
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(H  =  0.547-0.587),  with  only  Caqueta  slightly  lower  ( H  = 
0.481),  as  we  quoted  in  Results,  which  could  mean  that  this 
population  has  undergone  more  genetic  drift  than  the  other 
Colombian  populations  or  suffered  a  slightly  greater  founder 
effect  than  the  remainder  populations  analyzed.  Nevertheless, 
some  bottleneck  analyses  (the  theory  generated  by  Cornuet 
and  others40  and  Luikart  and  others41;  the  M  ratio  of  Garza  and 
Williamson42;  the  within  locus  kurtosis  (k)  test  and  the  interlo¬ 
cus  (g)  test  proposed  by  Reich  and  Goldstein43  and  Reich  and 
others44;  Kimmel  and  others’45  test;  Zhivotovsky  and  others’46 
test;  and  Beaumont’s  theory  with  the  MSVAR  software47) 
were  performed  (shown  elsewhere),  and  no  evidence  of  this 
phenomenon  was  observed.  Therefore,  this  loss  of  gene  diver¬ 
sity  was  more  likely  because  of  the  founder  effect  in  the  ori¬ 
gin  of  the  Colombian  populations  than  by  posterior,  or  recent, 
bottlenecks. 

A  comparison  of  gene  diversity  with  other  species  of  mosqui¬ 
toes  reveals  the  moderate  gene  diversity  found  in  Colombian 
An.  marajoara.  For  example,  Conn  and  others48  estimated  a 
value  of  H  =  0.834  for  An.  darlingi,  the  main  malaria  vector  in 
the  Americas.  Molina-Cruz  and  others49  determined  an  aver¬ 
age  H  =  0.78  for  An.  albimanus  in  Central  America  and  north¬ 
ern  South  America.  In  Africa,  An.  gambiae  also  yielded  higher 
gene  diversity  than  the  Colombian  populations  of  An.  mara¬ 
joara.  Lehmann  and  others50-51  estimated  values  for  different 
groups  of  populations  for  African  malaria  vectors,  and  the 
average  values  were,  respectively,  0.62  and  0.66.  In  an  iden¬ 
tical  sense,  other  studies  with  An.  gambiae  using  microsatel¬ 
lites  showed  values  ranging  from  0.60  to  0.85.52-53  The  same  is 
apparent  with  An.  maculatus  in  Thailand  ( H  =  0.83-0.92).54 

The  mosquitoes  sampled  in  the  current  study  were  obtained 
in  urban  and  rural  areas,  representing  diverse  ecological  con¬ 
ditions.  For  instance,  the  Norte  de  Santander  sample  was  col¬ 
lected  in  rice  fields,  whereas  the  Magdalena,  Antioquia,  and 
Meta  samples  came  from  more  urban  environments  and 
the  Caqueta  sample  was  obtained  in  an  Amazonian  habitat. 
In  other  studies,  populations  sampled  in  diverse  ecological 
situations  showed  significantly  different  gene  diversity  levels. 
Onyabe  and  Conn55  found  significant  heterozygosity  differ¬ 
ences  for  populations  of  An.  arabiensis  among  a  savannah  eco¬ 
system  (H  =  0.606-0.743),  a  transition  area  of  savannah-forest 
( H  =  0.642),  and  a  rain  forest  ( H  =  0.506).  However,  in  the  cur¬ 
rent  study,  no  heterozygosity  differences  were  found  among 
the  Colombian  populations  analyzed.  Thus,  natural  history 
and  demographic  events,  more  than  ecological  traits,  could  be 
the  main  facts  determining  the  gene  diversity  levels  found. 

With  a  few  exceptions,56  HWE  disequilibrium  by  homozy¬ 
gote  excess  is  typically  found  in  studies  of  malaria  vector  mos¬ 
quitoes  independently  of  the  molecular  markers  used.51-55-58  In 
the  present  study,  some  homozygote  excess  causes  could  be  dis¬ 
carded.  (1)  Null  alleles.  Although  it  is  never  totally  impossible 
to  discard  the  presence  of  null  alleles  using  PCR  procedures, 
the  deviations  from  HWE  for  each  one  of  the  markers  stud¬ 
ied  in  each  population  was  very  different.  If  null  alleles  were 
frequent,  we  could  expect  similar  homozygote  excess  trends 
in  those  loci  affected  by  this  phenomenon  in  different  popu¬ 
lations  if  the  genetic  heterogeneity  among  these  populations 
is  not  very  high  (as  it  is  in  the  current  case).  However,  with 
the  exceptions  of  M9-2  and  M6-8 ,  which  presented  homozy¬ 
gote  excess  in  all  the  populations  analyzed,  the  markers  stud¬ 
ied  showed  a  very  different  dynamic.  Therefore,  we  assumed  a 
weak  influence  of  null  alleles  on  the  homozygote  excess  found. 


(2)  Hitchhiking  and  synteny.  The  linkage  equilibrium  analysis 
only  revealed  one  case  of  significant  disequilibrium  between 
M9-2  and  M5-7  in  Meta  and  the  total  population.  However, 
none  of  other  loci  pairs  showed  any  evidence  of  significant  dis¬ 
equilibrium.  Possibly,  this  STRP  pair  could  be  placed  within, 
or  near,  a  chromosome  rearrangement  and  therefore,  is  sub¬ 
jected  to  adaptative  selection.59-60  Nevertheless,  with  the  excep¬ 
tion  of  the  aforementioned  case,  the  behavior  of  each  marker 
used  was  independent  from  the  other  loci.  This  absence  of  sig¬ 
nificant  gametic  disequilibrium  also  indicated  a  small  influ¬ 
ence  of  genetic  drift  in  the  populations  studied.  (3)  Endogamy. 
This  affects  the  whole  genome  identically.  Thus,  if  endogamy 
was  important  in  our  case,  we  would  expect  a  similar  homozy¬ 
gote  excess  for  all  the  markers  within  each  population  studied, 
which  we  did  not  detect.  Furthermore,  the  reproductive  biol¬ 
ogy  of  the  species  does  not  favor  endogamy.  The  ethological 
features,  such  as  mate  recognition,  multiple  oviposition  sites, 
blood-meal  location,  and  non-uniform  distribution  of  breed¬ 
ing  places,  make  the  mating  of  mosquitoes  of  the  same  breed 
unlikely.  Therefore,  endogamy  can  effectively  be  discarded  as 
a  possible  explanation.  (4)  Selection  in  favor  of  homozygotes. 
Selection  could  occur  in  one  locus  but  not  systematically  in 
the  majority  of  loci  studied  or  in  all  the  populations.  Thus,  the 
overall  homozygote  excess  detected  cannot  be  attributed  to 
selection.  In  summary,  microgeographic  Wahlund  effect  (sub¬ 
division  process  within  each  one  of  the  populations  studied, 
with  the  exception  of  Caqueta)  seems  to  be  the  most  probable 
cause  for  this  phenomenon.  The  existence  of  other  Anopheles 
taxa  in  our  sample,  which  could  produce  positive  homozy¬ 
gote  excess,  could  be  discarded  because,  as  we  previously 
commented,  all  the  exemplars  were  sequenced  for  ITS2  and 
white  genes  and  were  classified  such  as  An.  marajoara  with¬ 
out  any  doubt  (the  three  specimens  from  Magdalena  that 
did  not  belong  to  this  taxa  were  eliminated  from  the  analy¬ 
ses).  Thus,  the  picture  could  be  as  follows.  We  have  observed 
that  two  gene  pools  (Eastern  and  Western)  were  detected  in 
Colombia  partially  because  of  the  eastern  Andean  cordillera. 
Additionally,  within  each  population,  we  found  a  fraction  of 
individuals  with  typical  multi-locus  genotypes  of  other  popu¬ 
lations.  Only  a  small  fraction  of  these  migrants  were  hybrids 
between  these  Eastern  and  Western  gene  pools,  which  means 
that  the  majority  of  foreign  mosquitoes  in  a  given  population 
consisted  of  migrants  with  no  genetic  interchange  with  the 
original  mosquitoes.  This  produces  homozygote  excess  by  the 
Wahlund  effect,  even  in  very  restricted  local  areas.  This  has 
also  been  detected  in  other  malaria  vector  mosquitoes.51-61-63 
Significant  genetic  heterogeneity  (F  =  0.106)  between  two 
districts  of  Hermosillo  (Mexico)  was  discovered  for  Aedes 
aegypti ,M  Bosio  and  others65  showed  that  significant  genetic 
heterogeneity  was  present  in  Ae.  aegypti  populations  from 
Thailand  separated  by  25  km.  Similarly,  Herrera  and  others66 
found  very  restricted  gene  flow  among  Venezuelan  Ae.  aegypti 
populations  at  distances  of  15  km,  although,  in  these  areas, 
there  is  intensive  terrestrial  and  aerial  traffic. 

Although  two  gene  pools  were  detected  and  the  genetic  het¬ 
erogeneity  was  significant,  it  was  only  moderate  ( FST  =  0.069). 
Other  vector  insect  populations  have  presented  similar  genetic 
differentiation  levels  as  follows:  An.  gambiae  showed,  in  two 
Kenyan  populations  ( Jego  and  Asembo ) ,  a  value  of  F  =0.072- 
0.10,  although  the  genetic  differentiation  was  conspicuously 
lower  in  other  African  areas  separated  by  6,000  km  (FST  =  0.016).67 
Additionally,  Lehmann  and  others50  showed  that  the  genetic 
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heterogeneity  among  16  locations  in  10  African  countries  for 
An.  gambiae  ranged  from  FST  =  0-0.40,  with  a  mean  value 
around  FST  =  0.10.  Ruiz-Garcia  and  others6869  showed  similar 
genetic  heterogeneity  (F  =  0.067)  among  several  Colombian 
populations  of  Psoropliora  columbiae  in  a  geographical  area 
(Tolima,  Meta,  and  Cordoba  Departments)  that  overlapped 
with  our  study  area.  De  Souza  and  others70  showed  a  simi¬ 
lar  population  differentiation  in  Ae.  aegypti  from  Argentina 
( Fst  =  0.065)  as  well  as  An.  albimanus ,49  with  FST  values  from 
0.018  to  0.02  between  Central  American  populations  and  from 
0.057  to  0.082  between  South  American  populations,  similar  to 
those  shown  here.  Nevertheless,  other  species,  like  several  An. 
nuheztovari  populations  in  Colombia  and  Brazil,  presented 
substantially  greater  genetic  heterogeneity  ( FST  =  0.324), 71 
and  several  other  species  showed  an  underlying  lower  genetic 
heterogeneity,  such  as  Ae.  aegypti  from  Brazil  (Fs  =  0.018)72 
or  Ae.  albifasciatus  from  Argentina  (F  =  0.003), 73  than  the 
mosquito  populations  analyzed  here.  Therefore,  the  genetic 
heterogeneity  found  in  our  study  is  only  intermediate  with 
respect  to  values  found  in  other  mosquitoes  species.  It  is  pos¬ 
sible  that  human  activity  could  enhance  the  gene -flow  levels 
between  the  populations  of  An.  marajoara  studied,  which  has 
been  reported  for  other  mosquitoes  (Ravel  and  others74  for 
Ae.  aegypti  in  Mexico;  Vazeille  and  others75  for  Ae.  albopictus 
in  Madagascar;  Garcia-Franco  and  others76  for  Ae.  aegypti  in 
other  Mexican  populations;  Bosio  and  others65  for  Ae.  aegypti 
in  Thailand). 

The  migration  of  individuals  and  the  gene-flow  levels 
detected  were  elevated,  despite  the  fact  that  some  impor¬ 
tant  geographical  barriers  separated  the  populations  studied. 
Only  the  Eastern  Andean  cordillera  was  decisive  in  the  par¬ 
tial  differentiation  of  the  two  gene  pools  determined.  Some 
other  barriers  such  as  the  Cauca  and  Magdalena  rivers  in  the 
western  area  or  the  Macarena  Mountains  and  the  Tepuyes 
formations  in  the  eastern  region  studied  did  not  represent 
any  geographical  barrier  for  the  migration  of  An.  marajoara 
throughout  Colombia.  The  very  different  ecological  condi¬ 
tions  in  the  geographical  areas  analyzed  also  did  not  represent 
any  problem  for  the  dispersion  of  this  potential  malaria  vec¬ 
tor.  Malaria  has  recently  reemerged  as  a  public-health  burden 
in  Colombia.  From  18  to  24  million  people  are  living  in  areas 
with  high  risk  of  exposure  to  malaria,2  which  represents  about 
60%  of  the  total  Colombian  human  population.  In  2008  alone, 
around  100,000  cases  were  officially  reported  in  Colombia,77 
although  the  real  number  is  probably  higher.  Thus,  the  control 
of  the  malaria  vectors  such  as  An.  marajoara  is  very  important, 
because  it  shows  considerable  dispersal  and  a  flexible  adapta- 
tive  capacity  to  anthropogenic  environments.  An.  marajoara 
could  probably  be  playing  an  important  role  either  as  a  pri¬ 
mary  or  secondary  malaria  vector  in  some  Colombian  areas. 
The  occurrence  of  this  species  sympatrically  with  An.  darlingi 
adds  additional  complexity  to  the  study  of  the  dynamics  of 
malaria  transmission  in  Colombia.  From  a  population  genet¬ 
ics  point  of  view,  the  present  results  seem  to  suggest  that  the 
vector-control  programs  in  Colombia  do  not  affect  the  genetic 
structure  of  this  species  that  is  characterized  by  elevated  gene- 
flow  estimates  throughout  extensive  areas  of  this  country. 
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